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ABSTRACT 

We have performed cosmological simulations in a ACDM cosmology with and without radiative cooling, in 
order to study the effect of cooling on the cluster scaling laws. Our simulations consist of 4.1 million particles 
each of gas and dark matter within a box-size of 100 /;"' Mpc and the run with cooling is the largest of its 
kind to have been evolved to z = 0. Our cluster catalogues both consist of over 400 objects and are complete 
in mass down to ^ 10'"'/j~'Mo. We contrast the emission-weighted temperature-mass {Tg„-M) and bolometric 
luminosity-temperature (Lboi ~ 7e„) relations for the simulations at z = 0. We find that radiative cooling increases 
the temperature of intracluster gas and decreases its total luminosity, in agreement with the results of Pearce et al. 
Furthermore, the temperature dependence of these effects flattens the slope of the T^^i-M relation and steepens 
the slope of the Lboi - Tew relation. Inclusion of radiative cooling in the simulations is sufficient to reproduce the 
observed X-ray scaling relations without requiring excessive non-gravitational energy injection. 

Subject headings: galaxies: clusters: general, cosmology: theory 
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1. INTRODUCTION 

The mass of clusters of galaxies'* is dominated by dark mat- 
ter. The evolution of the dark matter halo population is now 
well-understood both theoretically (Lacey & Cole 1994) and 
numerically (Jenkins et al. 2001). The halos themselves are 
approximately self-similar and may be described in their inner 
regions by a one-parameter model (Navarro, Frenk & White 
1997) with a concentration parameter that is a slow function of 
mass. (Note, however, that there are significant deviations from 
this simple profile in the outer parts of clusters; see Thomas 
etal. 2001.) 

The intracluster medium (ICM) does not share the approx- 
imate self-similarity of the dark matter. This is expressed 
most clearly in the X-ray luminosity-temperature relation. For 
pure bremsstrahlung emission, the bolometric X-ray luminos- 
ity should scale with temperature as Lx oc (the inclusion of 
line emission flattens this relation), whereas observations in- 
dicate a much steeper temperature dependence, especially for 
low-mass systems (Edge & Stewart 1991; David et al. 1993; 
Ponmanetal. 1996; White, Jones & Forman 1997;Xue&Wu 
2000). At high temperatures part of this steepening is due to a 
central cooling flow, but removing the cooling flow component 
still does not reconcile the observations with the self-similar 
prediction (Allen & Fabian 1998; Markevitch 1998). 

The reason for the departure from self-similarity is that the 
gas is not as centrally-concentrated in clusters as the dark mat- 
ter, which is best physically expressed as an increase in entropy 
of the innermost gas (Evrard & Henry 1991; Kaiser 1991; 
Bower 1997). The most obvious explanation for this is that 
there has been some form of energy injection. The amount of 
energy required depends upon the density of the gas at the time 
when the heating occured. Lloyd-Davies, Ponman & Cannon 
(2000) argue for heating prior to cluster collapse and estimate 
a value of 0.3 keV per particle. Tozzi & Norman (2001) show 
that this can be lowered to 0.1 keV per particle by heating at 



the optimal time (when the gas is at its minimum density) but 
most other studies that consider heating within collapsed halos 
require much higher values of 1-3 keV per particle (Wu, Fabian 
& Nulsen 2000; Bower et al. 2000; Lowenstein 2000). 

There are two likely sources for any excess energy: stellar 
winds/supernovae and active galactic nuclei (AGN). It is known 
for CDM models that some feedback of energy into the inter- 
galactic medium must occur in order to prevent the "cooling 
catastrophe", in which the majority of the baryons in the Uni- 
verse cool and form stars at high redshift (White & Frenk 1991; 
Cole 1991;Blanchard, Valls-Gabaud&Mamon 1992). If flie 
heating efficiency is high, supernovae can inject an energy of 
order 0.3 keV per particle into the intergalactic medium, but 
they do not do so in an optimal way. Various authors, all of 
whom consider realistic, but different, models for the build-up 
of structure, conclude that supernovae are unable to provide the 
required excess entropy (Valageas & Silk 1999; Wu, Fabian 
& Nulsen 2000; Bower et al. 2000). Heating of the inter- 
galactic medium by quasars is also not without its problems 
as the heating must occur at just the right time in order not to 
overly suppress galaxy formation. Alternatively, the heating 
may arise from AGN buried within individual galaxies (Bower 
etal. 2000). 

Radiative cooling results in the removal of low-entropy gas 
in the cluster core and thus leads to an overall increase in tem- 
perature of the ICM (Thomas & Couchman 1992). The amount 
of cooling that takes place in a cooling flow after the final as- 
sembly of a cluster is insufficient to explain the observations 
(Bower et al. 2000), but a more realistic model in which cool- 
ing occurs at every stage of the collapse hierarchy can give en- 
tropy increases equivalent to an excess energy of 1-2 keV per 
particle (Wu, Fabian & Nulsen 2000) — in essence most of the 
cooling occurs in galaxy-sized halos before the formation of 
the cluster Bryan (2000) has developed a simple model in 
which low-entropy gas is removed from the cluster core and 
the surrounding gas is assumed not to have cooled at all. Al- 
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though this model has obvious deficiencies, it predicts the cor- 
rect luminosity-temperature relation from 0.5-lOkeV. 

Simulations of cluster formation including radiative energy 
loss have been carried out by Pearce et al. (2000). They showed 
that the gas is slightly heated and that the luminosity is greatly 
reduced (except in cooling flow clusters), in agreement with ex- 
pectations. However, the simulations were of limited resolution 
and covered only a small range in cluster mass. We are now un- 
dertaking a programme of simulations to extend these results 
over a wider mass-range and to contrast the properties of clus- 
ters in different cosmologies. In this letter we report results at 
z = from two simulations of a 100 hr^Mpc box in the ACDM 
cosmology, one with and one without radiative cooling. 

The simulations and cluster extraction method are described 
in Section 2 and the results are presented in Section 3. We sum- 
marize our conclusions in Section 4. 

2. METHOD 
2. 1 . The simulations 

We have carried out two simulations with 160^ particles 
each of gas and dark matter within cubical volumes of side 
100/i~'Mpc. The cosmological parameters were as foUows: 
density parameter, = 0.35; cosmological constant, Aq = 
A/3//2 = 0.65; Hubble parameter, h = i/o/lOOkms-iMpc"^ = 
0.71; baryon density parameter, ^Ibh^ = 0.019; power spec- 
trum shape parameter, T = 0.21; and a linearly-extrapolated 
root-mean-square dispersion of the density fluctuations on a 
scale 8 /!"'Mpc, erg = 0.90. With these parameters, the gas and 
dark matter particle masses are approximately 2.6 x 10^ and 
2.1 X IO^^H'^Mq, respectively. This mass is below the Stein- 
metz & White (1997) limit above which numerical heating 
dominates cooling. The runs were started at a redshift, z = 50 
and evolved to the present day, z = 0. The gravitational soft- 
ening was fixed at 50/j~^kpc in comoving co-ordinates until 
z = 1 and thereafter held constant at 25 hT^kpc in physical co- 
ordinates. This softening is sufficient to prevent two-body re- 
laxation (Thomas et al. 2001) and we have checked that the gas 
and dark-matter have similar specific energy profiles in clusters 
drawn from the non-radiative simulation. 

The only difference in the two runs was that one of them 
included radiative cooling. For this run, we used the cooling 
tables of Sutherland & Dopita (1993) and assumed a uniform 
but time- varying metallicity of Z = 0.3(1 /to)ZQ, where t/to is 
the age of the Universe in units of the current time. This time- 
varying metallicity is meant to crudely mimic the gradual en- 
richment of the ICM by stars, but more importantly it results 
in a global cooled gas fraction at the end of the simulation of 
approximately 20 per cent. This is the maximum value inferred 
from the observations of clusters (Balogh et al. 2001) which 
means that our simulation will represent the largest effect that 
cooling alone is likely to have on the intracluster medium. 

We use a parallel version of the HYDRA A^-body/SPH code as 
described by Couchman, Thomas & Pearce (1995) and Pearce 
& Couchman (1997) except that the SPH equations have been 
modified to use the pairwise artificial viscosity of Monaghan 
& Gingold (1983) — for a test of different SPH formalisms see 
Thacker et al. (2000). We decouple the hot and cold gas in the 
manner described by Pearce et al. (1999, see also Ritchie & 
Thomas 2001 ) to prevent artificial overcooling of hot gas onto 
the central cluster galaxies. Groups of 13 or more cold, dense 
gas particles (wifli S > 500 and T <12000K) are merged to- 
gether to form collisionless galaxy particles, which can only 



grow via the accretion of more gas. Not only does this save 
considerable computational effort, it also prevents small objects 
from being artifically disrupted within cluster potentials. 



2.2. The cluster catalogue 

Initially we identify clusters in our simulation by searching 
for groups of dark matter particles within an isodensity contour 
of 200, as described in Thomas et al. (1998). We work with a 
preliminary catalogue of all objects with more than 250 parti- 
cles, then retain only those which have a total mass within the 
virial radius exceeding Mum w 1.18 x \Q^^h~^MQ, correspond- 
ing to 500 particles of each species. The use of a small mass 
for the preliminary cluster selection ensures that our catalogue 
is complete. We have checked that using a different isodensity 
threshold, a different selection algorithm, or using gas particles 
instead of dark-matter particles to define the cluster, leads to an 
almost identical cluster catalogue — ^the only difference being 
the merger or otherwise of a small number of binary clusters. 

We define the centre of the cluster to be the position of the 
densest dark matter particle. Because the density parameter of 
the real universe is not known, we choose to average properties 
of the clusters within spheres that enclose an overdensity of 200 
relative to the critical density (whereas for this cosmology the 
virial radius corresponds to an overdensity relative to critical of 
about 110). Our final catalogues consist of 427 and 428 clusters 
in the radiative and non-radiative simulations, respectively. 



2.3. Cluster X-ray properties 
We calculate the bolometric luminosity of each cluster using 

L^i = ^j^MTi,Z), (1) 

where the subscript i denotes the sum over all gas particles 
within radius r2oo that have temperatures, 7] > lO^K, masses 
m, and densities, p,-; we assume a mean molecular mass /iwe = 
10~^'*g and an emissivity, A(7i,Z(f)), that is the same function 
used by HYDRA to calculate cooling rates, as discussed in Sec- 
tion 2.1. Note that, for temperatures below about lO^K, Une 
emission is important and the cooling rates are substantially 
higher than the contribution from bremsstrahlung alone. 

The emission-weighted temperature of each cluster is calcu- 
lated as 

_ J2imiPiHTi,Z)Ti 
2^,.OT,p,A(7;-,Z) 

Many of the smaller clusters in our catalogues have emission- 
weighted temperatures that are below 0.5 keV which means that 
the majority of their emission will emerge at energies that are 
below the X-ray bands. We do not attempt to calculate the emis- 
sion in any particular X-ray passband in this paper but instead 
quote bolometric luminosities and emission- weighted tempera- 
tures. 



3. RESULTS 
3.1. Temperature-mass relations 
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Fig. 1. — The emission-weighted temperature versus (total) mass within a 
density contrast of 200. Crosses and filled circles represent clusters from the 
simulations with and without radiative cooling, respectively. The various lines 
are described in the text. 

We plot the emission-weighted temperature-mass relation 
for both cluster samples in Fig. 1. Results from the simu- 
lations with and without cooling are illustrated using crosses 
and fiUed circles, respectively. The broken lines are power-law 
fits to the observational relation as determined by Homer et al. 
(1999), using mass estimates from galaxy velocity dispersions 
(dashed). X-ray temperature profiles (dash-dotted), the isother- 
mal /3-model (dotted) and the surface brightness deprojection 
method (dash- triple dotted). 

For the simulation with radiative cooling the relationship be- 
tween Tgyf and M is an approximate power law 



^7;w = i.9i 



/ M200 



V 1014/i-iMg 



0.58 



keV, 



and for the simulation without radiative cooUng 

0.67 



VT^ = 0.98 



M200 



keV. 



(3) 



(4) 



lOi4/j-iM0^ 

These should be compared with the virial relation (shown as the 
solid line on the figure) 

/ M200 



A;7;ir = 1.61 



V10i4/i-iMe 



keV. 



(5) 



(Note that no attempt has been made when making these fits 
to reproduce the observational selection effects. They should 
therefore be regarded as rough guides rather than precise pre- 
dictions.) The clusters from the non-cooling run mostly have 
emission-weighted temperatures that are much lower than the 
virial values. The reason for this is that their emission is dom- 
inated by high-density, low-temperature gas in the core of the 
cluster (this temperature drop in the core is simply due to the 
fact that the specific energy profile of the gas mimics that of 
the dynamically-dominant dark matter). Small fluctuations in 
the core properties of the clusters (which are not well-resolved 
by our simulations) lead to a large scatter in predicted temper- 
atures. It is important to note that these non-radiative simula- 
tions do not provide sensible predictions for the observed prop- 
erties of real clusters — the core gas has a short coohng time 
and would not in reality persist in the intracluster medium for 
the lifetime of the cluster. 



We note that our non-radiative Tew - M relation has a lower 
normalization than found by Thomas et al. (2001) for simu- 
lations in the rCDM cosmology because there it was assumed 
that the emission is purely bremsstrahlung, which underesti- 
mates cooling rates below lO^K and places less weight on the 
innermost region of clusters. 

In contrast to the non-radiative run, the clusters from the 
radiative simulation have emission-weighted temperatures that 
show less scatter and that exceed the virial values. The rea- 
son for this is that the core gas has cooled to low temperatures 
and been removed from the intracluster medium, leaving be- 
hind higher-entropy, higher-temperature gas — an effect that is 
more pronounced in lower-mass clusters. The emission is no 
longer dominated by the core and is well-resolved by our simu- 
lations. The clusters provide an adequate fit to the observational 
data given the large uncertainty in the latter as evidenced by the 
various broken lines in Fig. 1. 

3.2. Luminosity-temperature relations 
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Fig. 2. — The bolometric luminosity-temperature relations for our simula- 
tions. Clusters are represented by the same symbols used in Fig. 1 and broken 
lines represent the fits to observational data by Xue & Wu (2000). 

Fig. 2 illustrates X-ray luminosity-temperature relations for 
both simulations, again using crosses and filled circles for 
radiative and non-radiative clusters respectively. We have 
trimmed the original catalogues by selecting clusters only with 
temperatures above 0.35 and 0.75 keV for the non-radiative and 
radiative cooling catalogues respectively, to assure complete- 
ness in temperature. The broken lines illustrate best-fit power- 
law relations as determined by Xue & Wu (2000) for their col- 
lected group sample (dashed line), cluster sample (dot-dashed 
line) and combined sample (dotted line). (Note that the units for 
the normalisation of the relations in Table 2 of Xue & Wu 2000 
are misquoted and should be a factor of 10 larger, in agreement 
with their Table 1 and Figure 1 .) 

A power-law fit to results from the radiative simulation gives 

T 



Lbol = 9.0 X 10^ 



h ergs 



(6) 



Lbol = 1.2 X 10' 



44 I 



IkeVy 

and for the non-radiative simulation 

• 7; y 

J keV 7 

(As stated in section 3.1, these fits should not be taken as 
precise calibrations.) Again, both the slope and normalization 



-ergs 
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differ between the two simulations. The non-radiative simula- 
tion gives a slope close to 2, as predicted by self-similar scaling 
relations (the relation is slightly flatter due to the inclusion of 
line emission). However, the radiative simulation gives a signif- 
icantly steeper slope of 3.3. The normalization of the relation is 
around a factor of 10-100 lower in the radiative simulation than 
in the non-radiative simulation. Qualitatively, these differences 
are the same as were found by Pearce et al. (2000), for a smaller 
cluster sample. The high entropy gas that replaces cooled ma- 
terial in the radiative simulation is hotter and less dense than 
the gas in the non-radiative simulation. The change in density 
has the greater effect, since the X-ray emissivity is a slow func- 
tion of temperature but is proportional to the square of the gas 
density. The combination of the increase in temperature and 
the decrease in luminosity of the clusters causes the substantial 
shift in the Lboi - T^v/ relation. 

The results from the radiative simulation are in good agree- 
ment with the best-fit relations of Xue & Wu (2000). If any- 
thing, we predict luminosities that are are a factor of 2-3 too 
low, although there is stiU significant uncertainty in the obser- 
vational determinations. 

4. CONCLUSIONS 

In this letter, we have presented results from an ongo- 
ing programme to measure the evolution of X-ray cluster 
properties for a range of physical and cosmological mod- 
els. Specifically, we have presented the current-day (z = 0) 
emission- weighted temperature-mass {T^^i—M) and bolometric 
luminosity-temperature (Lboi - Tew) relations from two simula- 
tions of a ACDM cosmology, one with and one without radia- 
tive coohng. 

The Jew relation is significantly different in non-radiative 
and radiative simulations, with the latter in reasonable agree- 
ment with observational determinations. In the non-cooling 
simulation, the emission is dominated by cold, dense gas in the 
cores of the clusters; radiative cooling removes this gas from 
the intrac luster medium (converting it into stars) and replaces it 
with higher-entropy, hotter material. This effect is more preva- 
lent in lower-mass systems and so flattens the temperature-mass 
relation. 

The Z,boi - T'ew relation is also significantly different between 



the two simulations. The high entropy material in the radia- 
tive simulation is less dense than the material it replaces and 
causes the X-ray luminosity of clusters to decrease by around 
a factor of 100 at Tew = 1 keV. The slope of the relation in the 
non-radiative simulation is 1.9, similar to self-similar predic- 
tions. However, the slope of the radiative relation is signifi- 
cantly steeper, 3.3, again due to the differential effect of cool- 
ing with halo temperature. The radiative simulation is in much 
closer agreement with the observations, both in the slope and 
normalization of the relation. 

In this paper, we have used bolometric luminosities and 
emission-weighted temperatures. For the low-temperature clus- 
ters found in the non-radiative simulations, these will differ 
significantly from properties measured in any particular X-ray 
band. However, we wish to emphasize that to attempt to correct 
for this is misleading as the use of such simulations is wrong 
in principle — ^the low-entropy gas in the cores of these clusters 
has a short cooling time and will not be present in real systems. 

The global fraction of cooled gas (and stars) in the radiative 
simulation, 20 per cent, is higher than suggested by observa- 
tions of the A'-band galaxy luminosity function (Balogh et al. 
2001). More importantly, the fraction of cooled gas within the 
clusters varies between about one-third and two-thirds with de- 
creasing mass. While these values are not convincingly ruled 
out by observations, most people would also regard these as 
high values — ^in which case our simulation can be treated as an 
upper bound on the effect of radiative cooling. 

In this paper, we have deliberately ignored the effect of non- 
gravitational heating upon the gas. In reahty, we know that 
there must be heating associated with star-formation and metal- 
enrichment of the intracluster medium. This will raise the en- 
tropy of the gas and reduce the amount of cooling that is re- 
quired to match the observations. However, we do not regard 
the case for significant heating by AGN or very efficient super- 
novae feedback as proven. 

The simulations described in this paper were carried out on 
the Cray-T3E at the Edinburgh Parallel Computing Centre as 
part of the Virgo Consortium investigations of cosmological 
structure formation. OM is supported by a DPST Scholarship 
from the Thai government; PAT is a PPARC Lecturer FeUow. 
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